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Summary. — 

Highly variable 7-ray pulses and X-ray flares in GRB light curves can result from 
external shocks rather than central engine activity under the assumption that the 
GRB blast-wave shell does not spread. Acceleration of cosmic rays to > 10 20 eV 
energies can take place in the external shocks of GRBs. Escape of hadronic energy 
in the form of UHECRs leads to a rapidly decelerating GRB blast wave, which may 
account for the rapid X-ray declines observed in Swift GRBs. 

PACS 95.85.Ry - 98.70.Rz. 



1. Introduction 

GRB light curves measured with Swift consist of a BAT light curve in the 15 - 150 
keV range followed, after slewing within as 100 s, by a detailed 0.3 - 10 keV XRT X-ray 
light curve [T]. This information supplements our knowledge of the highly variable hard 
X-ray and 7-ray light curves measured from many GRBs with BATSE and other GRB 
detectors. About one-half of Swift GRBs show X-ray flares or short timescale structure, 
sometimes hours or later after the onset of the GRB. Approximately 30% of the Swift 
GRBs display rapid X-ray declines, and an additional w 30% display features unlike 
simple blast wave model predictions [2]. 

We make three points in this paper: 

1. Highly variable light curves can be produced by an external shock under the as- 
sumption that the GRB blast wave does not spread, or spreads much more slowly 
than assumed from gas-dynamic or relativistic hydrodynamic models that do not 
take into account magnetic effects in GRB blast waves. If this assumption is valid, 
then it is wrong to conclude that highly variable 7-ray emissions, X-ray flares with 
At/t -C 1, or late time X-ray flares require delayed central engine activity or col- 
liding shells. 

2. External shocks in GRB blast waves can accelerate cosmic ray protons and ions 
to £ 10 20 eV, making GRBs a logical candidate to accelerate the highest energy 
cosmic rays. 



© Societa Italiana di Fisica 



1 



2 



CHARLES D. DERMER 



3. Escape of ultra-high energy cosmic rays (UHECRs) takes place from an external 
shock formed by an expanding GRB blast wave on time scales of a few hundred 
seconds for the observer. Blast-wave deceleration due to the loss of the internal 
hadronic energy is proposed [3] to be the cause of X-ray declines in GRB light 
curves observed with Swift. 

2. X-ray flares and 7-ray pulses from external shocks 

We have performed a detailed analysis of the interaction between a GRB blast- wave 
shell and an external stationary cloud [4]. The analysis is performed under the assump- 
tion that the cloud width A c / <C x, where x is the distance of the cloud from the GRB 
explosion. The interaction is divided into three phases: (1) a collision phase with both 
a forward and reverse shock; (2) a penetration phase where either the reverse shock has 
crossed the shell while the forward shock continues to cross the cloud, or vice versa; and 
(3) an expansion phase, where both shocks have crossed the cloud and shell, and the 
shocked fluid expands. The shell width is written as 

(1) A(x) = A + 77 -| , 

1 

and the proper number density of the relativistic shell is given by 

E 



(2) n(x) 



4:TTx 2 TQm p c 2 A(a 



where Tq is the coasting Lorentz factor of the GRB blast wave, and Eq is the apparent 
isotropic energy release. 

Short timescale flaring requires (a) a strong forward shock, which from the relativistic 
shock jump conditions [5] imply a maximum cloud density given by 

(3) n cl < * 



167rx 2 rQm p c 2 A(a;) ' 

and (b) significant blast-wave deceleration to provide efficient energy extraction, which 
occurs in clouds with thick columns [BJ, that is, with densities 

(4) n cl > * 



These two conditions translate into the requirement that 
(5) A d > 4r 2 A(x) 



in order to produce short timescale variability. The short timescale variabilty condition 
[6J for quasi-spherical clouds is 
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Using eq. fl} for the shell width, eqs. ([5]) and ([6]) imply the requirement that 
(7) rj < 1/4T 

in order to produce rapid variability from an external shock. Hence the production of 
7-ray pulses and X-ray flares from external shocks depends on whether the GRB blast- 
wave width spreads in the coasting phase according to eq. (fT]), with rj Rj 1, as is generally 
argued. In the gas-dynamical study of [7], inhomogeneities in the GRB fireball produce 
a spread in particle velocities of order \v — c\/c ~ r^ 2 ' so that A(x) ~ x/Tq when 
x ^ TqAo- This dependence is also obtained in a hydrodynamical analysis [5]. 

Two points can be made about these relations. First, the spread in A considered 
for a spherical fireball is averaged over all directions. As the fireball expands and be- 
comes transparent, the variation in fluid motions or gas particle directions over a small 
solid angle ~ 1/Tq of the full sky becomes substantially less. Second, the particles 
within a magnetized blast-wave shell will expand and adiabatically cool so that the fluid 
will spread with thermal speed v t h = PthC The comoving width of the blast wave is 
r A + PthcAt' ps r A + Pth x/Tq, so that the spreading radius x spr = T^A n /(3 t h- Adia- 
batic expansion of nonrelativistic particles can produce a very cold shell with (3q ;$ 10~ 3 , 
leading to very small shell widths. 

The requirement on the thinness of A(x) does not apply to the adiabatic self-similar 
phase, where the width is necessarily ~ x/Tq, as implied by the relativistic shock hy- 
drodynamic equations [5]. Even in this case, however, A -C x/Tq if the blast wave is 
highly radiative [9]. Under the assumption of a strong forward shock and small clouds 
in the vicinity of a GRB, highly variable GRB light curves are formed with reasonable 
efficiency ( 10%) to transform blast wave energy into 7 rays [6} ITU]. 



3. Cosmic ray acceleration in GRB blast waves 

The maximum particle energy for a cosmic ray proton accelerated by an external shock 
in a GRB blast wave is derived. Consider a GRB blast wave with apparent isotropic 
energy release Eq — 10 54 i^54 ergs, (initial) coasting Lorentz factor Tq = 300r3oo, and 
external medium density uq — IOO712 cm~ 3 . The comoving blast wave volume for the 
assumed spherically symmertric explosion, after reaching distance x from the center of 
the explosion, is 

(8) V = 4ttx 2 A', 

where the shell width A' = x/12T (the factor 1/1 2r is the product of the geometrical fac- 
tor 1/3 and the factor 1/4T from the continuity equations of relativistic hydrodynamics; 
r is the evolving GRB blast wave Lorentz factor). 

The Hillas condition |TT] for maximum particle energy E' max is that the particle 
Larmor radius is less than the size scale of the system; E max in the stationary frame 
(primes refer to the comoving frame) is given by 

(q \ J _ E 'max _ E rnax w 

1 ' ~ eB> TeB> ' 

The largest particle energy is reached at the deceleration radius x = x& when T = r , 
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where the deceleration radius 

(10) x d = (— ^2 f' 3 = 2.6 x io 16 (-^-)V3 cm . 

Hence E max = ZeB'x d /l2 . 

The mean magnetic field B' in the GRB blast wave is assigned in terms of a magnetic 
field parameter es that gives the magnetic field energy density in terms of the energy 
density of the downstream shocked fluid, so 



(11) B' = (32 7 rnoeBTOpC 2 ) 1/ V r ( r - 1) = 0.4(e B n ) 1/2 r 1200Ve^r 3 oo Gauss 



Thus 

(12) E rnax S 8 x lQ^Zn^e^Tl^Eli 3 eV 

[T2l |T5] . so that external shocks of GRBs can accelerate particles to ultra- high and, 
indeed, super-GZK energies. Implicit in this result is that acceleration occurs within the 
GRB blast wave through, for example, second-order Fermi acceleration [13]. Acceleration 
to ultra-high energy through first-order relativistic shock acceleration requires a highly 
magnetized surrounding medium [14] . 



4. Rapid X-ray declines from UHECR escape 

If UHECRs are accelerated by GRB blast waves, then blast-wave dynamics will be 
affected by the loss of internal energy when the UHECRs escape. This effect is proposed 
to explain the rapid X-ray declines in the Swift GRB light curves [3]. Photohadronic 
processes become important when the threshold condition e'7' Si m n /m e cr ~ 400, where 
e = hv/m e c 2 is the dimensionless photon energy, m p c 2 j is the proton energy, and 7 is the 
proton Lorentz factor. For protons interacting with photons at the peak photon energy 
e p k — 2Te' pk /(l + z) of the vF v spectrum, 

™ p ^, 3 x 10 16 (r/300) 2 
(!3) E pk 77- — ; eV. 

(1 + z)e p k 

The comoving timescale for a proton to lose a significant fraction of its energy through 
photohadronic processes is given by ^(7), where t'^(^) ~ {K^a^jle'n'^e')]^ K^a^ 
= 70 fib is the product of the photohadronic cross section and inelasticity, and the co- 
moving energy density of photons with energy s» e' is u' e , = m e cre n' ph (e'). 

The relation between the measured vFv flux f e and internal energy density is u' e , = 
d\f t / (ca; 2 r 2 ), where = 10 28 c?28 cm is the luminosity distance of the GRB. For protons 
interacting with photons with energy e' fc , we therefore find that the comoving time 
required for a proton with energy E p k (as measured by an observer outside the blast 
wave) to lose a significant fraction of its energy through photohadronic processes is 

(u , ., (F m e c 2 x 2 r 2 e; fc x 2 6 (r/300)(l + z) epfc 

where x = 10 16 £i6 cm and f e k = 10 _6 /_6 ergs cm -2 s _1 is the vF v flux measured at 
e p k ; the relation between E p k and e p k is given by eq. (fT3"]) . 




Observer time t(s) Observer time t(s) 

Fig. 1. - Rates and inverse timescales as a function of observer time for 10 20 eV cosmic ray 
protons as measured by a stationary external observer. Left and right panels are results for 
parameter sets 1 and 2, respectively. 



The dependence of the terms x(t), f €pk (t), T(t), and e p k(t) on observer time in eq. (fTi|) 
can be analytically expressed for the external shock model in terms of the GRB blast 
wave properties Eq, Tq, environmental parameters, e.g., uq, and microphysical blast 
wave parameters cb and e e [3]. This can also be done for other important timescales, for 
example, the (available) comoving time t' ava since the start of the GRB explosion, the 
comoving acceleration time t' acc = C, acc mpC 2 ^' / eBc, written as a factor Q acc 3> 1 times the 
Larmor timescale |15j . the escape timescale t' esc in the Bohm diffusion approximation, 
and the proton synchrotron energy loss timescale t' syn . 

Fig. 1 shows the rates (or the inverse of the timescales) for 10 20 eV protons in the 
case of an adiabatic blast wave that decelerates in a uniform surrounding medium. The 
left-hand panel of Fig. 1 uses the parameter set 

z = 1, Tq = 300, £ 5 4 = 1, n = 1000 cm" 3 , e e = 0.3, e B = 0.3 , 

and the right-hand panel uses the parameter set 

z = 1, T = 150, E54, = 10, n = 1000 cm" 3 , e e = 0.1, e B = 0.3 . 

The characteristic deceleration timescale in the left and right cases, given by td — 9.6(1 + 
z)(E 5 4 : /n2Tl 00 ) 1 / 3 s, is w 9 s and w 120 s, respectively. 

For these parameters, it takes a few hundred seconds to accelerate protons to en- 
ergies w 10 20 eV, at which time photohadronic losses and escape start to be impor- 
tant. Photohadronic losses inject electrons and photons into the GRB blast wave. The 
electromagnetic cascade emission, in addition to hyperrelativistic electron synchrotron 
radiation from neutron escape followed by subsequent photohadronic interactions [16j . 
makes a delayed anomalous 7-ray emission component as observed in some GRBs [171118) . 
Ultra-high energy neutrino secondaries are produced by the photohadronic processes. De- 
tection of high-energy neutrinos from GRBs would confirm the importance of hadronic 
processes in GRB blast waves. The ultra-high energy neutrons and escaping protons 
form the UHECRs with energies £ 10 20 eV. 
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The GRB blast wave rapidly loses internal energy due to the photohadronic processes 
and particle escape. The blast wave will then rapidly decelerate, producing a rapidly 
decaying X-ray flux. As argued in more detail elsewhere [3], the rapidly decaying fluxes 
in Swift GRBs are signatures of UHECR acceleration by GRBs. If this scenario is 
correct, GLAST will detect anomalous 7-ray components, particularly in those GRBs 
that undergo rapid X-ray declines in their X-ray light curves. 

* * * 
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